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Abstract: Astronomical light echoes, the time-dependent light scattered by dust in the vicinity of 
varying objects, have been recognized for over a century. Initially, their utility was thought to be 
confined to mapping out the three-dimensional distribution of interstellar dust. Recently, the discovery 
of spectroscopically-useful light echoes around centuries-old supernovae in the Milky Way and the Large 
Magellanic Cloud has opened up new scientific opportunities to exploit light echoes. 
In this review, we describe the history of light echoes in the local Universe and cover the many new 
developments in both the observation of light echoes and the interpretation of the light scattered from 
them. Among other benefits, we highlight our new ability to spectroscopically classify outbursting 
objects, to view them from multiple perspectives, to obtain a spectroscopic time series of the outburst, 
and to establish accurate distances to the source event. We also describe the broader range of variable 
objects whose properties may be better understood from light echo observations. Finally, we discuss 
the prospects of new light echo techniques not yet realized in practice. 

Keywords: ISM: reflection nebulae — ISM: supernova remnants — stars: individual (RS Puppis, 
S CrA, R CrA, Nova Persei 1901, V838 Mon, 77 Car) — supernovae: individual (SN 1987A, Cas A, 
Tycho) — ISM: individual (SN 1987A, Cas A, Tycho, SNR 0509-675, SNR 0519-69.0, SNR N103B) - 
novae — stars: distances 



1 Introduction 

Modern astrophysics has long benefited from the addi- 
tional information conveyed by variable stars. Some 
of the most physically interesting sources are those 
which emit pulses of light during eruptive or disruptive 
events which traditionally have resulted in brief tem- 
poral windows during which critical physical informa- 
tion about rare, transitional phases can be acquired. 
In recent years, the ability to extend those temporal 
windows through detecting and analyzing light echoes 
from historical events has been realized and a num- 
ber of eruptions without contemporaneous photometry 
and spectroscopy have become observable with mod- 
ern instrumentation. 

"Light echoes" (hereafter, LEs) are simply faint re- 
flections of the light from a bright astrophysical source 
off material - in these cases, circumstellar or interstel- 
lar dust. The sizes and irregular surfaces of typical 
dust grains results in incident light being scattered 
over a wide range of angles although forward-scattering 
is typically far stronger than back-scattering. The 
scattered light received by an observer on Earth has 
taken a longer path than any light which may have 
been received directly from the outbursting source, and 
so it arrives later. This time-delay depends on the 
geometry betw een the dust location and the source 
event (Section 3.1 1. As we shall describe presently, 



the extreme brightness of supernovae (SNe) have al- 
lowed detection of their LEs many centuries after the 
observed event. In the intervening time, we have de- 
veloped physics and our instrumentation capabilities 
to a point where we can gain tremendous new insights 



on the nature of these rare events - a situation which 
would be impossible without the temporal delay intro- 
duced by LEs. 

This review reports and summarizes the wealth of 
new information which has become available as a result 
of studying LEs from varying astrophysical sources, fo- 
cusing on scattered light echoes in the local Universe. 
It is now over a century since t he first LEs were discov - 
ered around Nova Persei 190 1 (jRitchey|1901a|b||T90"2"l), 
and re cognized as such by Kapteynp 1902 1; Perrine 



(1903). Initially, LEs were mainly used to constrain 



the 3D-position of the scattering dust and its prop- 
erties like grain-size distribution, density, and compo- 
sition. The impact of unresolved light echoes on the 
spectra and light curves of SNe was recognized as a 
complication in the analysis of extra-galactic SNe. We 
now recognize that currently observationally-difEcult 
or impossible problems of asymmetry estimation, spec- 
tral typing, lightcurve reconstruction, and more pre- 
cise distances yield to LE observations. From the wealth 
of applications listed above, we emphasize in this re- 
view the more recent techniques which use LEs as a 
means to directly observe the outburst light of an event 
by using a non-direct line-of-sight to the source. The 
full potential for the use of LE techniques has not yet 
been realized and we identify key exploration oppor- 
tunities for the near future. 



2 Finding Light Echoes 

The first LEs were disco vered around Nova Persei 1901 
(Ritchey 1901a|b 1902[ ), and it was interpreted as such 
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by Kapteyn (19021; Perrine (1903 k Since then, LEs 



have been seen associated with a wide variety of ob- 
jects: the Galactic Nova Sagittarii 1936 ( |Swope|1940| 



the e ruptive variable V838 M onocerotis (|Bond et al 



2003 b, the Cepheid RS Puppis (WesterlundTl96T 
19721, the T Tauri star S CrA (Ortiz et al.' 



and the Herbig Ae/Be star R CrA (Ortiz et al 



Havlcn 



2010 



2010 



sients could be found, and subsequent targeted searches 
found LEs of Tycho's SN flRest et al.||2007| |2008b |, 
Cas A (Rest et al.||2007| |2008b| jKrause et al.||2008a[ ), 
and i) Carinae (Rest et al.| 2012 1 . A recent search for 
LEs of historic SNe around 4 recent SNe in M83 using 



Echoes have also been observed from extragalactic SNe, 
with SN 1987A being the most famous case (Crotts| 



polarization imaging was not successful (Romaniello 
It should be noted 



et al 



(2005 



20051 



that 



Krause et al 



1988 S untzeffet al. 1988b), but a lso including SNe 1980 K 
(Sugerman et al.||2012| |, 19 91T (Schmidt et al. 



Sparks et al.||1999| ), 1993 J (Sugerman fc Crotts 



1994 



2002 



iu et al.||2003|), 1995E (Quinn et al.||2006|> , 1998bu 
( Garnavich et al.|2001[|Cappellaro et al.|2001[ ), 2002hh 
( Welch et al.||2007| |Otsuka et al.||2012[), 2003gd (Sug-| 
|erman|2005| |Van Dyk et al.|2006| |Otsuka et al.|2012[) 



2004ct (Otsuk a et al.|2012[ >, 200 6X (Wang et al. 2008 
Crotts fc Yourdon|2008| 72006 bc (Gallagher et al.]2011 
Otsuka et al.||2012[), 2006gy (Miller et 



discovered "IR echoes" from the Cas A SNR, in 
which the dust absorbs the outburst light, is warmed, 
and re-radiates light at longer infrared (IR) wavelengths 
This is different to the usual LE phenomenon described 
here, where the light is simply scattered by dust, there- 
fore preserving the spectral characteristics of the source 
event. 

LEs are extended objects, often with very faint sur- 
face brightnesses. Therefore it is necessary to apply 



difference imaging (Tomaney fc Crotts||1996| |Alard fc 



al.||2010[ ) 

2007it (Andrews et al.|2011| >. One factor that all these 
objects have in common was that the LEs were de- 



Lupton 



sky bac 



1998| |Alard||2000[ ) to separate them from the 



(ground in order to be able to identify them 
and measure their surface brightness. However, only 
the relative fluxes between two epochs is revealed in a 
difference image. Therefore obtaining absolute fluxes 
for individual epochs has traditionally relied on a single 
template image that is free of LEs, which often must 
be constructed by a complicated and usually subjective 
process of hand-selecting suitable images (e.g., [&uger-| 
man et al.||2005b| ). |Newman fc Rest| pOOo) presented 
a solution to th i s issue by applying the NN2 method 
observed for' SN 199T (Schmidt et al.||1994[ ), a SN la of|Barris et~al~1 p005| to extract the relative fluxes 



tected serendipitously in the presence of the still lu- 
minous stars or transients. These light echoes may 
influence the spectra as well as the light curves of the 
sources when unresolved, in particular for events like 
type II SNe that are likely located in dust-rich envi- 
ronments (e . g., |Schaefer|1987b| |Chugai|1992| |Di Carlo| 



|et al. 



2002; Otsuka et al. 



20121, but have also been 



that is not necessa rily expected to be in such a dust- 
rich environment. Roscherr & Schaefer (20001 finds 



that light echoes cannot explain the extremely long 
decline of SN Iln like SN 1988Z and SN 1997ab. 

The suggestion that historical SNe might be stud- 
ied by their scattered LEs was first made by |Zwicky| 
( 19401. Simple scaling arg uments (Shklovskii 1964 
van den Bergh 1965bj |1966|) based on the visibility of 



Nova Persei predicted LEs from SNe as old as a few 
hundred to a thousand years can be detected, espe- 
cially if the illuminated dust has regions of high density 



(> 10" 



The few dedicated surveys in the last century for 
LEs from historic SNe (van den Bergh 1965a b 



Boffi et al.||1999|l and novae (van den Bergh 



1966 



1977 



Schaefer|1988[ ) have been unsuccessful. However, these 
surveys did not use digital image subtraction tech 
niques ( Tomaney fc Crotts 1996 Alard & Lupton 1998 



of LEs across a range of epochs directly from a series 
of difference images. This method treats all images 
the same and makes maximal use of the observational 
data. The efficacy of this technique was demons trated 
by applying it to the LEs around SN 1987A (New-| 
|man fc Rest|200"6| ). The most common source of false 
LE candidates is scattered/reflected light from bright 
stars falling in the focal plane beyond the edges of 
the detectors. Small pointing errors between image 
epochs can produce difference features whose shapes 
and surface brightnesses can mimic LEs. These fairly 
common optically-induced difference features make it 
difficult to identify faint LEs with software and most 
LE searches rely at least in part on visual inspection 
of the difference images. With the advent of the next 
generation of wide-field, time domain surveys like Pan- 
STARRS (Kaiser et alT1[2010l, PTF (Rau et al.|2009l ), 



|Alard|[2000| " to remove the dense stellar and galactic 
backgrounds. Even the bright ech oes near SN 1987A 
(Crotts|l988| 



Skymapper (Keller et al. 20071, and ultimately LSST 



(Ivezic et al. 2008[ ), visual inspection will no longer be 



Suntzeff et al 



1988b I at V ~ 21.3 mag arcsec 



a viable option and software solutions for identifying 
-strue light echoes will need to be developed. 



are hard to detect relative to the dense stellar back- 



ground of the Large Magellanic Cloud (LMC). Maslov 



( 2000[ ) suggested that using a wide-field polarization 3 Light Echo Primer 



imager to detect LEs of historic SNe. 

The situation changed with the advent of CCDs 
and telescopes with large field-of- views, which enabled 
the astronomical community to entertain the first wide- 
field time-domain surveys with sufficient depth. The 
first LEs of ancient events were found serendipitously 
by |Rest et~aT[ (2005b|) as part of the Super MACHO 
survey ( |Rest et al.||2005a| : They found LEs associ- 
ated with 400-900 year-old SNRs. This result demon- 
strated that the LEs of historical SNe and other tran- 



There are many excellent reviews of LEs covering the 
various aspects of LE science. The groundwork of 



LE geometry was laid by Couderc (19391, and the 
surface brightness of LEs was derived by |Chevalier] 
(1986| ). After SN 1987A started to show its beauti- 
ful set of LEs, more derivations o f the surface bright- 
ness of single-scattered LEs (e.g., |Schacfcr 1987a| [Xu| 
et al.|1994 Sugerman|2003| |Patat|2005|l and mu ltiple- 
scattered LEs ( Patat||2005| |Patat et ah] [2006 ) were 
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done, as well as its impact on observed SN spectra 
and light curves discuss ed (e.g.,|Schaefer|1987b| |Patat| 
|2005||Patat et al.|2006"| . |Sugerman| ( |2003| investigated 
the range of circumstances which would produced ob- 
servable light echoes for different classes of transients. 

Besides the scattered light echoes which preserve 
the SED of the source event, there are other types of 
time-delay which are often referred to as echoes. One 
of them is the so-called "IR echo" , where dust absorbs 
the outburst light, is warmed, and re-radiates light at 
longer infrared (IR) wavelengths. IR echoes have been 
observed around SNe (e.g.jBode fc Evans|1980b| Dwek 



1983 



Krause et al.|2005||Kotak et al. 2009 ) and novae 
[agT jBode fc Evans||1980a| |1985| |Gehrz||1988"l ). "Re- 
combination" echoes are another form of echoes where 
the initial light gets absorbed and then re-radiated at 
a different wavelength (e. g., |Panagia et al.|1991 Gould 
1994||Gould fc Uza|1998| . The ionization light echo of 



a quasar in "Hanny's Voorwerp" was discovered by the 
Galaxy Zoo project ( Lintott et al.| 2009 Rampadarath 
et al.|20T0 i. All of the above are different to the usual 



LE phenomenon described here, where the light is sim- 
ply scattered by dust, therefore preserving the spectral 
characteristics of the source event. Because IR and 
recombination echoes only share the geometry with 
scattered light echoes, we will focus on scattered light 
echoes, and only refer to the other echoes when they 
are relevant to a particular object with scattered light 
echoes. 

3.1 Geometry 



Figure [T| (fr om |Sugerman (2003), adapted from Xu 
|et al.| ( |1994[ )) illustrates the geometry of the LE phe- 
nomenon. Light originating from the SN event source, 
O, is scattered by dust located a distance z in front 
of the SN and a projected distance p perpendicular 
to the line of sight, and is redirected towards the ob- 
server. The event source and observer are separated 
by a dist ance D. In app lication z « D, and the LE 
equation ( Couderc|1939 1 



2rf 



ct 
~2 



(1) 



can be derived, where t is the time since the explosion 
was originally observed and c is the speed of light. 
Then the distance r from the scattering dust to the 
source event is r 2 = p 2 + z 2 , and the projected distance 
on the sky is p — (D — z) tan 7, where 7 is the angular 
separation between the source event and the scattering 
dust. 

Measuring p through imaging, we can determine 
the exact three-dimensional location of the scattering 
dust by simply knowing the time since the outburst oc- 
curred and the distance to the object. In the case that 
the distance to the object is not well known, the above 
geometry still provides relative distances between dis- 
tinct scattering dust locations. Through LE imaging 
alone, the three-dimensional dust structure in front 
of an outburst event can be mapped in great detail. 
While the dust structure can be measured by observ- 
ing multiple scattering dust locations, it is most thor- 
oughly mapped through monitoring the apparent mo- 




Figure 1 : Figure from Sugermanl (12003b, which was 
adapted from Xu et al. '{ 1994 ). Note that distances 



are not to scale and D is much larger than z. 



tion of a given LE system, since the apparent motion 
is heavily dependent on the scattering dust inclination 
(see Section [3~3]>. 



3.2 Surface Brightness 

The surface brightness of LEs has bee n derived in dif- 
ferent variants by various autho rs (e.g.,|Chevalier|1986 
Schaefer|1987al |Xu et al |1994| |Sugerman|2003| |Patat 



it 



2005 1 . Here we follow a derivation of the surface bright- 
ness flux fsB by Sugerman (2003J and define 



f SB (X,t,4,) = F(X)n{T)(^^S{X,e). 



(2) 



This equation has four principal components: (1) F(X) 
is the lightcurve-weighted integrated event flux of the 
varying source event, (2) n(r) is the dust density as 
a function of position r, (3) a wavelength-independent 
spatial component, where Ap and Az are the width 
of the LE on the sky and the depth of the scattering 
dust sheet along the line of sight, respectively, and (4) 
the wavelength-dependent integrated scattering func- 
tion, S(X, 8), which we describe in more detail in Sec- 
tion ETO 

Note that the event flux F(X) is the flux differ- 
ence of the source with respect to its quiescent state. 
In all the cases where the source event is a one-time 
transient (e.g. SNe) or where the event is several mag- 
nitudes brighter than the magnitude in the quiescent 
state (e.g. V 838 Mon), this differentiation is not im- 
portant since the quiescent flux is either zero or close 
to zero. However, for other objects like Cepheids this 
is important since what we can detect is the appear- 
ing and disappearing scattered flux which is due to the 
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brightness variation (i.e. the difference in flux between 
its maximum and minimum brightness) and not the to- 
tal brightness. For example, we can find scattered light 
from a non-varying source at a given position, but it is 
constant with time (i.e. there is no apparent motion of 
the scattered light) and is commonly called a reflection 
nebula. 

It is important to note that this derivation assumes 
a photon is only scattered once (single-scattering ap- 
proximation). However, for very dense dust environ- 
ments multiple scatter ing can become important ( Chev i- 
lier|1986| |Patat|2005[ ). The effect that unresolved sin- 
gle and even more so multiple scattered LEs have on 
observed colors an d spec tra of SNe is discussed in de- 
tail inlPatat et al.|(|2006|. 



3.2.1 Dust Scattering 

The integrated scattering function S(A, 9) describes 
the wavelength-dependent effect of the scattering off 
dust grains, and is therefore of great importance when 
calculating the surface brightness of LEs. Here we de- 
scribe how S can be calculated for our Galaxy and the 
Magellanic Clouds. 

Since different types of dust grains have different 
integrated scattering functions, they must be added 
together to get the total integrated scattering function 
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Figure 2: The upper panel shows the integrated 
scattering function S(X, 0) for the three dust mod- 
els "MWG", "LMC av g", an d "SMC bar", which 
are defined in Section 3.2.1| The bottom panel 
shows S(X, 6) normalized by S(8OOuA,0) for the 
three models and various scattering angles. 



where X denotes the dust grain type. Here, X can 
represent silicon dust grains, carbonaceous dust grains 
with a neutral Polycyclic Aromatic Hydrocarbon (PAH) 
component, or carbonaceous dust grains with an ion- 
ized PAH component ( Weingartner & Drainc 2001 1. 

The integrated scattering function Sx for a dust 
grain of type X is 



Sx(X,6) = J Qsc,x{K a)a g §x{9,\, a)fx{a)da 



(4) 



where Qx (A, a) is the dust grain scattering efficiency 
and fx (a) is the dust grain density distribution. The 
radius of the individual dust grains is given by a, with a 
grain cross-section a g — iva 2 . The Henyey-Greenstein 
phase function, $x, is given below, with gx(^,a) de- 
noting the degree of forward sca ttering for a given 
grain (Henyey & Greenstein||1941 1. 



$x(A,0) 



Si (A, 



(l + g% (A, a) - 2g x (A, a) cos 0) 3 / 2 



(5) 



Values for Qx and (A, a) can be derived using tables 
provided by B. T. Drains 1 ( Draine fc Lee|l984||Laor fc 



Draine|199 3 Weing artner fc Draine|200lj|Li fc Draine 



2001) 



We use the "MWG", "LMC avg", and "SMC bar" 



models defined by Weingartner & Drainc (2001), adopt- 
ing values for their model parameter be of 5.6 x 10~ 5 , 
2 x 10 -5 , and 0, respectively (Table 1 and 3, Weingart- 



|ner fc DTa inc 2001). The models consist of a mixture 



of carbonaceous grains and amorphous silicate grains. 
The dust grain size distribution is 



/(«) 



1 dn 



gr 
da ' 



(6) 



1 http: / /www. astro. princeton.edu/~draine/dust/dust. did 



where n gr (a) is the number density of grains smaller 
than size a and tlh is the number density of H nuclei 
in both atoms and molecules. Note that small car- 
bonaceous grains (a < 10"Vm) are PAH-like while 
large carbonaceous grains (a > 10 _3 /xm) are graphite- 
like ( |Li fc DTa inc 2001). The size distributions for 
carbonaceous dust (fd{a) — C'ionf(a) and f C n(a) — 
(1 — C'ion)f(a)) can then be calculated using equations 
2, 4, and 6 of Weingartner & Drainc (2001). Assum- 
ing PA H/graphitic grains t o be 50% neutral and 50% 
ionized ( Li fc DrainepOOl I, don =0.5. Equations 5 
and 6 are used for amorphouse silicate dust, f s . 

In the upper panel of Figure [2] we show the in- 
tegrated scattering function S(X, 9) as a function of 
wavelength for the Galaxy and Magellanic Cloud mod- 
els, adopting a scattering angle 9 = 90°. In the 4000A to 
8000A region shown, the Milky Way Galaxy dust scat- 
tering is almost twice as efficient in the blue. For a 
detailed comparison of the scattering efficiencies, we 
show in the bottom panel of Figure [5] S(X, 9) normal- 
ized by S(9>0QQA,9) for the three models and various 
scattering angles. Note that the difference in blueward 
scattering efficiency is much larger between 9 = 20° 
and 60° than between 60° and 150°. The differences 
between the models are small, but significant, and do 
not seem to vary significantly with scattering angle, 
htmln Figure [3] we show ^(A, 0) in log-scale normalized 
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Figure 3: S(X, 9) normalized by S{\,90°) for 
the three dust models "MWG", "LMC avg " and 
"SMC bar", which are defined in Section 3.2.1 
Large differences in models and wavelengths are 
only apparent for scattering angles 9 < 60°. 



by S(X, 90°) for the different models and wavelengths. 
For scattering angles 9 > 60°, the scattering efficiency 
is quite similar for all models and wavelengths. For 
6 < 60°, however, the efficiency of blueward scattering 
is more prominent, with large differences between dif- 
ferent models and different wavelengths clearly visible. 



3.3 Apparent motion 

One of the most easily measurable properties of LEs 
is their (often superluminal) apparent proper motion 
on the sky. The main influences on the apparent mo- 
tion are the angular separation between the LE and 
source event, the time since the first light of the source 
reached earth, and the distance to the source event (9, 
t, and D, respectively). This opens up the possibility 
that if either one of the age or distance to the event is 
known, the other parameter can be estimated since the 
angular separation can be measured very accurately. 
However, the apparent motion also depends on the in- 
clination of the scattering dust filament with respect 
to the line of sight. This fact is often overlooked or 
ignored, and therefore the derivation of age/distance 
of the event assumes implicitly or explicitly a certain 
dust inclination. This can lead to an underestimate 
of the systematic uncertainty and subsequently to a 
wrong scientific conclusion. A tale of caution is the 



analysis by Krause et al. ( 2005 I of the 1R echoes of 



Cas A. Their main scientific conclusion is that most if 
not all of these IR echoes are caused by a recent X- 
ray flare of the Cas A SNR based on their apparent 
motions. However, the analysis was flawed because it 
did not account for the fact that the apparent motion 
strongly depends on the inclination of the scattering 
dust filament ( |Dwek fc Aren dt 2008; Rest et al., in 
prep.). 

Figure [4] illustrates how different dust sheet incli- 
nations produce different apparent LE motions. We 
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Figure 4: The ellipses of equal arrival time defined 
by Equation [l] are shown as the solid and dashed 
black line for t — 300 and t = 301 ly, respectively. 
The dust sheets A, B, C, and D have inclinations 
of a = Q, 45, 90, and 135 degrees, respectively. 
LEs scattering off dust sheets A, B, C, and D 
have apparent motions of 1.0c, 0.5c, 0.0c, and oo, 
respectively. Here, c is the speed of light. 



define the inclination angle, a, as the angle of the dust 
sheet with respect to the p axis, where positive an- 
gles go from the positive p axis towards the negative 
z axis. The ellipses of equal arrival time defined by 
Equation [I] are shown as the solid and dashed black 
line for t = 300 and t = 301 ly, respectively. The dust 
sheets A, B, C, and D have inclinations of a = 0, 45, 
90, and 135 degrees, respectively, crossing the point 
(z,p) = (0,300) ly. LEs scattering off dust sheets A, 
B, C, and D have apparent motions of 1.0c, 0.5c, 0.0c, 
and oo , respectively. Here, c is the speed of light. One 
can see that for dust inclinations in Figure [4] where 
—45 < a < 0, the apparent motion of the LE is su- 
perluminal. Theoretically, for each apparent motion 
there exists one unique dust filament inclination in a 
180 degree range that produces exactly that apparent 
motion. For a LE at a given position, there exists 
degeneracy between the time t since explosion, the 
dust inclination, and the apparent motion. In order 
to determine one of these parameters, the other two 
parameters must be known. However, the inclination 
of the dust filament is seldom known. A solution is to 
marginalize over the dust inclination, showing that the 
apparent motion associated with a dust filament with 
an inclination perpendicular to the LE ellipsoid is a 
very good approximation of the expectation value (see 
Rest et al., in prep., for a detailed derivation). Note 
that this marginalization assumes that the inclination 
of the dust is random, which might not be true due to 
detection biases. As an example, a common miscon- 
ception is that for z = 0.0, i.e. the dust is in the plane 
of the sky, the apparent motion of LEs is c, implicitly 
assuming that the dust inclination is a = 0.0°. How- 
ever, the true expectation value is v = 0.5c associated 
with a = 45°. As we will see below, the inclination is 
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an important physical property of the scattering dust 
that not only affects the apparent motion of the LE, 
but the observed flux profile and spectrum of the LE 
as well. 

3.4 Geometric Distance 

Under suitable circumstances, LEs provide a unique 
opportunity to determine the geometric distance to its 
source object. For a given light echo at an angular 
distance to the source event with known event time, 
there is still a degeneracy between distance D (source 
to observer) and the distance z (scattering dust to the 
source event along the line of sight). This degeneracy 
gets broken if z can be estimated with some additional 
information or assumptions. 

One of the easiest assumptions to make is z — 0, 
or that the apparent motion of the light echo is the 
speed of light. However, if these assumptions are not 
based on any additional evidence, the systematics can 
be very large and consequently can lead to incorrect 
distances. Examples for this are discussed for V838 
Mon (Section IP) . 

If the geometrical shape of the scattering dust can 
be constrained, the degeneracy can also be broken. 
This has been done in the case of SN I987A where 



Panagia et al. ( 1991 1 determines the distance to the 



LMC to a few percent by assuming that its circum- 
stellar ring is circular (Note, however, that they use 
gaseous emission lines, and not scattered light echoes, 
which light travel delays follow similar equations) . |Gould| 
(19941); iGould & Uzal dl998b show that the distance 



changes only on the percent level if the ring is moder- 
ately elliptical. 

Another complication is added if the source of the 
light echoes has recurrent maxima. In that case it 
is difficult to determine which light echo corresponds 
to which maximum. If the wrong association is done, 
catastrophic errors can be the result, and therefore ex- 
tra care needs to be applied when determining these 
associations. We discuss two cases in Section [4.11 and 

EP1 

One of the most accurate and robust ways of de- 
termining the distance with light echoes is using the 
polarization of the scattered light. The scattering po- 
larizes the light, where the polarization of this scat- 
tered light will vary strongly with the scattering angle 
between the source and the observer (Draine 20031. 



For normal scattering cross sections, the angle of max- 
imum polarization, ip po i, is close to 90° and depends 
only slightly on dust properties like type and grain size 
distribution. If the time since the source event, t, is 
known and ^) vo i = 90° , the distance to the source event 
can be calculated with ct = Dtan(7), where 7 is the 
angular separation of the maximum polarization to the 
source event. In the event that ip po i is not exactly 90°, 
but is known, the method can still be used and the 
corrected distance is 



D' = Dx 



cos(V> po i) + 1 
sin(Vvoi) 



(7) 



distances. 

Since the maximum polarization angle ip po i depends 
on the type of scattering dust (usually unknown), sys- 
tematic uncertainties on the order of 5%-10% are in- 
troduced into these distance measurements (see Fig. 
5 of[t>rame 2003). For example, at a wavelength of 
6165 A, ippoi can range from 90° to 97° for SMC 
bar dust and Galactic dust, respectively. Therefore by 
determining ip po i at different wavelengths we can ver- 
ify if the polarization is consistent with the standard 
Galactic Rv = 3.1 dust model. If the dust is indeed 
consistent with Rv = 3.1, then the theoretical predic- 
tion for ippoi is well constrained and the systematics 
will be lowered to a few percent. 

The most important ingredient for this method 
is to have LEs with scattering angles spread around 
90°, which corresponds to having LEs with distances 
spread around 2 = 0. The best example where this 
method has been applied is V838 Mon, as discussed in 
Section [4. 3| which has an extensive light echo system 
spanning a wide range in z. However, not all objects 
with LEs have such extensive light echo systems, but 
Tycho flRest et al.||2008b[ ), Ca s A ([Rest et al.||2008b| ), 
and most recently n Carinae (Rest ct al. 2 012[ ) might 
have a suitable number of LEs that cover the required 
range of angular separations from the source event for 
this method to work. 



3.5 Light Echo Profile 

The flux profile of a LE is the slice through the LE 
along the p axis pointing toward the source event. This 
profile is actually the lightcurve of the source event 
stretched or compressed depending on the inclination, 
a, of the scattering dust filament, and convolved with 
the effects of the dust width, cr^, and the seeing (Rest 
et al.|2011b l 



As in Rest et al. ( 2011b 1, we illustrate this stretch- 



We refer the reader to Sparks (1994 1996) who first 



investigated how this method can be used to determine 



ing and compressing effect using a realistic LE scenario 
similar to that observed in the Cas A echo system. 
We consider an echo originating from an event in the 
Galaxy 300 years ago at a distance of 10000 light-years. 
We assume the outburst event was similar to the ob- 
served outburst of SN 1993 J. The double-peaked na- 
ture of SN 1993J's lightcurve makes for a particularly 
good example for illustrating the observational effects 
of the scattering dust and the PSF size. The left pan- 
els in Figure [5] show scattering dust filaments (brown- 
shaded area) of different widths and inclinations. The 
blue-shaded area indicates a 130 day long event simi- 
lar to SN 1993J with a fast rise, a peak after 30 days, 
and a slow decline. Note that the beginning of the 
outburst is observed at the largest distance away from 
the source on the sky, not the closest. The scatter- 
ing dust filaments are located at (z, p) = (0, 300) ly 
with inclinations of a = (0,0,0,45)° and widths of 
cr d = (0.008, 0.03, 0.2, 0.008) ly from top to bottom, 
respectively. We project the coordinate p onto the sky 
and show on the y-axis the angular difference in p to 
the peak in arcsec. 

The middle and right panels of Figure [5] show the 
LE flux profiles (black lines) for the different dust con- 
figurations and PSF sizes. For an infinitely thin dust 
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sheet and an infinitely small PSF, the flux profile is 
just the projected light curve, as indicated by the red 
line. For very thin dust filaments (cr<j = 0.008 ly) 
with a = and HST-like PSF size, the observed LE 
profile still shows the signature of the double peak 
in the lightcurve (upper middle panel). For ground- 
based seeing conditions (upper right panel), however, 
the two peaks are smeared out and indistinguishable. 
Similarly, if the dust width increases, more and more 
of the original lightcurve shape is lost (ad = 0.03 ly, 
middle panel in the 2nd row), until not much of the 
original shape is left (era = 0.2 ly, middle panel in the 
3rd row). Note that the effects of dust width and PSF 
size are somewhat degenerate. 

The effect of the dust filament inclination is illus- 
trated in the bottom panel, where the LE flux profile 
is squeezed by a factor of two compared to the pro- 
file shown in the top row due to the inclination of the 
dust filament of a = 45°. Fortunately, the dust fila- 
ment inclination is straightforward to determine from 
the apparent motion if the date of and distance to the 
source event is known ( Rest et al.|2011b l. A more de- 
tailed discussion of the interplay between dust width, 
dust in clination, and seeing can be found in |Rest et al.| 
(2011b I. Such accurate modeling of the dust is essen- 



tial when looking for small differences in spectral fea- 
tures, e.g., when comp aring spect r a of L Es from the 
same object as done in Rest et al. (2011a I. 



3.6 Light Echo Spectroscopy 

While imaging LEs can reveal the nature of the ma- 
terial around an outburst event, and possibly the dis- 
tance to the event (as in Section 3.4|, it is difficult to 
further our understanding of the outburst itself through 
imaging. Spectroscopy of LEs, however, allows the 
outburst to be studied in detail after it has already 
faded. The first time a spectrum of a LE was taken 
was in 1902 of Nova Persei in a 35 hour effort that can 
only be called heroic. It confirmed that the nebulous 
moving features seen around Nova Persei were indeed 



its echo es. Gouiffes et al. ( 1988 ) and Suntzeff et al. 
( |1988b[ ) obtained spectra of the inner and outer ring 
of SN 1987A's light echoes, and compared to averaged 
spectra from SN 1987A. They found the LE spectra 
were most similar to those of the SN near maximum 



light. Serendipitously, Schmidt et al. ( 19941 found that 



750 days after maximum, the spectrum of SN 1991T 
showed its light echo spectrum reflected from fore- 
ground dust superimposed on the late-time nebular 
spectrum. Since 2008 with the work of |Rest et al.| 
(2008a I, studies using LE spectroscopy have shown it 
to be a powerful tool to spectroscopically type ancient 
and historic SNe. 

Initial studies assumed that an observed LE spec- 
trum represented a lightcurve-weighted integration of 



the SN spectra at individu al epochs (e.g., Rest et al. 
2008a Krause et al.|2008a l. In most cases, this is ad- 



equate for an approximate spectral identification and 
classific ation of the source event. However, Rest et al. 
( |2011b| > showed that an observed LE spectrum repre- 
sents an integration weighted with an effective lightcurve. 
The effective lightcurve can be constructed by applying 



a window function to the originally observed lightcurve. 
Modeling both astrophysical (dust inclination, scatter- 
ing, and reddening) and observational (seeing and slit 
width) effects for a given LE are required to determine 
the correct window function and therefore correctly in- 
terpret the observed (integrated) LE spectrum. This is 
illustrated in Figure[5] If a slit is placed on the LE pro- 
file, in the cases where the width of the dust filament 
is not very large, the window function is narrow and so 
the slit probes only a limited part of the source event 
epochs (1st and second row of panels from the top). 
Only for very thick dust filaments (3rd row of panels 
from the top) is the window function sufficiently wide 
such that the spectrum probes all or nearly all epochs 
of the source event. 

An important part of determining the correct win- 
dow function is determining if late-time spectral fea- 
tures are included in the observed LE spectrum. Naively, 
one might expect these late-time features are not sig- 
nificant since the surface brightness of the source event 
already significantly declined. However, at late times 
when the spectrum is nebular, the observed flux is typ- 
ically concentrated into a few strong lines. We demon- 
strate the importance of this late-time contribution 



using LE spectroscopy from SN 1987A (Rest et al. 
2011b| ). The SN 1987A system is the ideal testbed 
for LE spectroscopy, where the high-flux LEs can be 
paired with the complete spectral and photometric his- 
tory of the SN as it was originally observed. 

Figure [6] shows a modified and upda ted version of 
Figure 18 and 19 in Rest et al. (2011b I, which com- 
pares a SN 1987 A LE spectrum to its model. The 



left panel shows the lightcurve of SN 1987A (|Hamuy 
et al.||T988l |Suntzeff et"al~1|1988a| |Hamuy fc Suntzeff 



1990) and the effective lightcurve resulting from mul- 



tiplying the lightcurve with the appropriate window 
function. The window function is determined from 
the best-fitting modeling of the LE profile, taking dust 
width and inclination, seeing and slit width into ac- 
count. The necessity of this modeling is clear when 
comparing model spectra (colored lines) to the ob- 
served LE spectrum (black line) in the middle and 
right panels of Figure [6] The model spectra shown 
here are integrations of the originally observed spectra 



of SN 1987A ||Menzies et al. 


1987; 


Catchpole et al. 


1987 


|1988| |Whitelock et al. 


1988 


Catchpole et al. 


1989 


Whitelock et al.|1989 Phillips et al.| 1988 19901 



line) lightcurves shown in the left panel. Here, mod- 
eling the correct window function limits the contri- 
bution of late-time spectra into the final integrated 
model spectrum. In the case of SN 1987A, the late- 
time spectra show strong nebular emission flux in Ha. 
The Ha line is shown in the right panel of Figure |6j 
where it is fit well with the effective lightcurve model. 
Incorrectly modeling the observed LE spectrum using 
the full lightcurve of SN 1987A results in an emission 
component to the Ha line that is far stronger than ac- 
tually observed in the LE spectrum. It is clear from 
Figure [6] that detailed consideration of the scattering 
dust and observational effects are critical in correctly 
interpreting LE spectra. Examples of LE spectroscopy 
are discussed in Section|4.5| 14.61 14.71 and 14.81 
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Figure 5: The blue-shaded area indicates a 130 day long event similar to SN 1993J with a fast rise, 
a peak after 30 days, and a slow decline. The scattering dust filament (brown shaded area) is located 
at ( z , p) = (0, 300) ly with dust widths <j d of 0.008 ly, 0.03 ly, 0.2 ly, and 0.008 ly, from top to bottom, 
respectively. All dust filaments have an inclination of a = 0°, with the exception of the bottom panel 
which has an inclination of a — 45° . The corresponding LE flux profiles are shown in the middle and right 
panels for PSF sizes of 0.05" and 1.0", respectively. For an infinitely thin dust sheet and an infinitely 
small PSF, the flux profile is just the projected light curve, as indicated by the red line. The phase of this 
projected light curve is shown on the right y-axis. 
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Figure 6: Analysis for SN 1987A LE spectrum taken using the Gemini Multi-Object Spectrograph on 
Gemini-South. Left: Original lightcurve of SN 1987A compared with the effective lightcurve associated 
with the observed LE. Middle: Observed LE spectrum of SN 1987A (black line) compared with the modeled 
spectra created by integrating the original outburst spectra of SN 1987A with the two lightcurves shown in 
the left panel. Right: Closeup of Ha line, showing that only the modeled effective lightcurve can correctly 
account for the relative strength of the Ha line in the observed spectrum. 




Figure 7: 3D illustration of Cas A LEs. North is toward the positive- y axis (up), east is toward the 
negative- X axis (left), and the positive- z axis points toward the observer with the origin at the Cas A 



SNR (r ed filled circle). The four brown circles indicate the scattering dust of LEs discovered by Rest et al. 
(2008b I. The black lines show the path of the light scattering from the LE-producing dust concentrations. 
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3.7 3D Spectroscopy 

LEs offer the exciting opportunity to obtain spectra of 
the original source event from different lines of sight 
(LoSs). Fig. [7] illustrates how LEs from Cas A's SN, 
which are scattered by different dust structures, probe 
the SN from different directions. Observing these arcs 
is equivalent to observing different hemispheres of the 
photosphere, providing direct observations of potential 
asymmetries in the source event. |Smith et al.| ( |2003a| ) 
first applied this technique, observing the n Carinae 
central star from different directions using spectra of 
the reflection nebula. 

The signatures of asymmetries in spectra of the 
same events are much more subtle than when doing 
spectral classification by comparing different events. 
Therefore it is essential that both astrophysical (dust 
inclination, scattering, and reddening) and observa- 
tional (seeing and slit width) effects are taken into 
account accurately, as described in the previous Sec- 
tion [3^6] and in more detail in |Rest et al.| pbllb[ ). An 
example of 3D-spectroscopy is discussed in Section |4~7l 



3.8 Spectroscopic Time Series 

One as-yet unrealized opportunity provided by SN LEs 
is the ability to obtain a spectral time series of the 
source event of LEs using very thin scattering dust 
filaments, as suggested by Rest et al] ( |2011b| ). As dis- 
cussed in Section |3.5[ the LE profile is simply the spa- 
tially projected light curve of the source event, con- 
volved with the finite thickness of the scattering dust 
filament and the finite size of the PSF. The spatial ex- 
tent of the projected lightcurve is determined by the 
inclination of the dust filament, in combination with 
the duration of the event. For short events like SNe, 
and typical dust configurations in our Galaxy, the spa- 
tial extent is on the order of arcseconds (see Figure [5| . 
This projected lightcurve is then smeared due to the 
thickness of the dust filament, in combination with 
the PSF size. Whether it is possible to temporally re- 
solve the SN spectra depends on the combination of 
the above parameters 

Figure [5] illustrates this: if the dust filament is very 
thin, for exa mple as observed for Cas A LEs (Rest 
et al.|2011b| ), the resulting window function can have 
a width of ~10 days. In the example of SN 1993J, a 
SN lib similar to Cas A, it would allow us to observe 
the spectrum of the shock breakout with an instru- 
ment with HST-like PSF-size (top and bottom middle 
panel). An enormous advantage of using LEs for ob- 
serving the shock break-out is that there is no urgency 
to obtain observations. Every pixel is a data point, 
and pixels before the edge of the LE are equivalent to 
closely spaced non-detections in a SN search. With 
this method, getting a spectrum at early epochs is as 
easy (or difficult) as getting a spectrum at late epochs, 
which is in stark contrast to SN searches for which it 
is difficult to get spectra at early epochs due to the 
lag time between when the image is taken, the SN is 
discovered, and follow-up is triggered. 

Spectra at early epochs of a SN or other transients 
are especially valuable scientifically, since they contain 



many signatures of the original explosion and also the 
progenitor. For example, with the shock break-out, the 
duration, luminosity, and SED depends on a handful 
of parameters such as the presence of a stell ar wind 
and the ejecta mass (Matzner & McKee||1999 1, but is 



most dependent on the progenitor radius ( Calzavara & 
Matzner|2004[). The peak brightness of the subsequent 



fading is dependent on the ejecta mass, kinetic energy, 
and progenitor radius, while the timescale for the fad- 
ing is proportional to both the radius of the photo- 
sphere and its temperature ( jWaxman et al.|2~0 07). The 
ability to temporally resolve the SN spectra allows one 
to follow the evolution of the explosion, e.g., the evo- 
lution of the velocity gradients of spectral features and 



the abundances of elements in the ejecta JRest et al. 



2011b| 



A particularly promising case for this technique is 
if the duration of the source event is very long, since 
the spatial extent of the projected lightcurve is pro- 
portional to its duration. For example, since the Great 
Eruption of r\ Carinae spanned more than a decade, it 
is now possible to obtain LE spectra of it from different 
epochs (Rest et al. 20121, which are only marginally 



affected by dust width and PSF-size of the observa- 
tions. 



4 Light Echo Case Studies 

In recent years, new LEs have been discovered for a 
variety of events, and new techniques have been used 
to utilize them. In this section we discuss some of 
the most prominent LE systems and note the achieve- 
ments, difficulties, and shortcomings of their analysis 
(see also Table [T]). 

4.1 S CrA & R CrA 



In a recent paper, |Ortiz et al.| ( |2010| ) presented the re- 
sults of their imaging of LEs around two young stellar 
objects, the T Tauri star S CrA, and the Herbig Ae/Be 
star R CrA. These classes of variable stars are known 
for irregular (non-repeating) variability. The authors 
draw attention to the discovery of the reflection nebu- 



losity and its variation by Hubble ( 1922 1 and Lightfoot 



( 1989 1 although the early work did not interpret their 



findings as LEs. 

Th e observational material available to lOrtiz et all 
( |2010[ ) for its analysis of S CrA consisted of I-band 
photometry from the All Sky Automated Survey ( |Po-| 
|jmanski||1997| spanning approximately 150 days and 
nine images spanning 93 days. They adopted a model 
of a spherically-symmetric shell of dust centered on 
the star and then associated the four brightness peaks 
from the ASAS time-series with the observed LE fea- 
tures from their imaging. Based on their assumptions 
they arrived at a distance to S CrA of 128 ± 16 pc, 
a distance between the star and dust shell of 10 4 AU. 



Ortiz et al. (20101 speculated that the dust shell could 



be similar in nature to the Oort cloud surrounding our 
Sun. 

The major uncertainty in distance estimates us- 
ing LEs is the scattering geometry which is unknown 
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without polarization measurements. Distances derived 
with polarimetry are also subject to systematic errors 
as we shall describe in Section |4.2| and [4. 3| even when 
the dust distribution is described correctly. The anal- 
ysis of S CrA LEs has several additional flaws. First, 
given four brightness peaks in four months and nine 
images spanning three months, it is very likely that 
there is at least one repeating LE - that is, a LE ap- 
pearing twice caused by two different peaks scattering 
off the same dust filament. Second, the brightness dif- 
ferences between the peaks do not seem to be large 
enough to be able to explain this. Third, there is a 
large analysis degeneracy between the distance to the 
star, the times since peak brightness(es), and the dis- 
tance from the dust to the star. This degeneracy has 
not yet been fully explored. Finally, the time series it- 
self is noisy and the identification of brightness peaks 
is not unambiguous. While it may be possible to de- 
termine a useful geometric distance to S CrA using LE 
methods and polarimetry, we conclude that the exist- 
ing distance estimate is flawed. 

4.2 RS Pup 

RS Puppis is one of the longest-period and consequently 
most luminous classical Cepheids known in our Galaxy. 
As such it is one of the Galactic counterparts of the 
bright Cepheids used to determine distances to the 
most distant galaxies amenable to distance estimation 
with this standard candle. Since Cepheids are believed 
to be a high-reliability distance indicator, determining 
an accurate, independent "geometric" distance to RS 
Puppis, which at a distance of approximately 2.0 kpc 
cannot be determined by trigonometric parallax with 
current instrumentation, would be highly desirable. 

The reflection neb ula of RS Puppis was discovered 
by |Westerlund ( |1961 1, who remarked on the possibility 
of LEs due to the variable nature of the star. IHavlenl 



( 1972 1 used a series of photographic images to demon- 



strate that features within the nebula show light vari- 
ations at the Cepheid's pulsation period and argued 
that these "echoes" could be used to derive a geometric 
distance, provided that the phase angle of the scattered 
light was known. More recently, [Kervella et al. ( 2008 1 



obtained an observing season's worth of high-angular 
resolution CCD images and used phase lag measure- 
ments of the RS Puppis LEs to estimate a geometric 
distance of 1992 ± 28 pc, which if correct, would be 
by far the most accurate distance to a Cepheid. ft 
would also be a distance unaffected by uncertainty in 
reddening. 

However, Bond fe Sp arks (2009J showed that the 
analysis by Kervella et al.| ( |2008[ ) has serious flaws. 
They pointe d out that the implicit assumption by |Kerv^ lla 
et al. ( 2008[) that the scattering dust lies, on average, 



in the plane of the sky is not valid since the scattering 
efficiency strongly favors forward scattering and thus 
significantly biases the results. In addition, |Bond fc| 
Sparks] 12009b showed that the IKervella et al.l (12008b 



results have low statistical significance and argued that 
with polarization measurements, as demonstrated with 
V838 Mon by |Sparks et al.| ( |2008| , the true geometric 
distance might yet be determined with high precision 



using this technique. 



4.3 V838 Mon 



Some of the most spectacular LEs ( |Henden et al.|20"02 



|Bond et al.||2003| l appeared when V838 Monocerotis 
went through an explosion in 2002 ( |Brown et al.|2002 1 , 
an outburst so energetic that it became one of the 
brightest stars in the Local Group for a few weeks in 
2002 at Mv = —10, gaining 9 magnitudes in the V 
band from its typical quiescent brightness. It ejected 
so much debris that the material is still not optically 
thin and is now thought of as a member of a new class 
of variables, the so-called intermediate-luminosity red 
transients. 

The echoes around this star were subsequently used 
to determine its distance. The first attempts using LE 
expansion measurements greatly underestimated the 
distance to V838 Mon to be < 1 kpc ( |Munari et al.| 



|2002 Kimeswenger et al. 20021, because they were 



based on the assumption that the LE apparent mo- 
tion is at the speed of light. More realistic estimates 
were first done using polarization measurements from 



HST imaging, setting a lower limit of < 6 kpc (Bond 



et al 



2003 1. A detailed construction of the 3D map 
of the scattering dust by |Tylenda] ( |2004[ ) provided a 
similar limit of < 5 kpc, and also indicated that the 
scattering dust is of interstellar and not circumstellar 
nature. Assuming that the dust is in a thin sheet, 
Crause et al.| ( |2005[ ) finds a distance of 8.9 ±1.6 kpc. 



The most a ccurate distance det ermination to date is 
provided by |Sparks et al.| ( |2008 l, whose detailed anal- 
ysis of the LE polarization and expansion speed yield a 
distance of 6.1 ± 0.6 kpc, in very good agreement with 
other methods (e.g., |Afsar fc Bond|2007| >. 

The LEs also show a remarkable feature caused by 
a double helix in the scattering dust filaments, which 



point almost radially towards V838 Mon (Carlqvist 



2005). A possible mechanism for the formation of such 



a complex structure is the twisting of outflowing ma- 
terial by a magnetic field, as suggested by |Carlqvist| 
(20051. 



V838 Mon also shows IR echoes (Banerjcc et al 



2006), which are spatially coincident with the scattered 



Mq) makes 



LEs. The large mass of the dust (> 0.2 
it unlikely that the dust is circumstellar and from a 
previous mass-loss episode, supporting the hypothesis 
that the dust causing both the scattered and IR echoes 
is interstellar ( [Banerjee et al.| |2006). 



4.4 SN 1987A 

Shortly after SN 1987A exploded, spectacular LEs scat- 
tered off interstellar dust at distances of ^100 and 
~400 pc (|Crotts|[l988l [Suntzeff et al.||1988b| |Gouiffes 
|et al.|1988||Couch et al.|l990p . Later, LE s scattered off 
circumstellar dust wer e also discovered (Crotts et al. 
(1989] |Bond et al.|1989[ ). These LEs have revealed valu- 
able information about the circumstellar environment 
around SN 1987A ( jSugerman et al.|[2005a|b[), the in 



terstellar medium (ISM) of the LMC (Xu et al.||1994 



Xu & Crotts|1999b, and about the asymmetry of 
WT987A itself (Sinnott et al., in prep). 
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4.4.1 SN 1987A light echoes scattered off 
circumstellar dust 



Crotts & Kunkel ( 1991 I uses the light echoes within 
10" of SN 1987A to determine the density and 3D 
structure of the scattering dust, constraining the mass 
loss of SN 1987A's progenitor. They found that the 
scattering dust is consistent with being from mass lost 
at a constant rate from a red supergiant atmosphere. 
Later observations revealed a double-lobed nebula, with 
a waist that is nearly coincident with the elliptical 
circumstellar ring seen in atomic recombination lines 
flCrotts et al]p95t. 



Sugerman et al. (2005a I uses LEs detected within 
30" of SN 1987A to construct the most detailed 3D- 
model of the circumstellar environment. They found 
that a richly structured bipolar nebula surrounds SN 1987A. 
with an outer, double-lobed "peanut" extending 28 lt- 
yr along the poles, which is believed to be the contact 
discontinuity between the red supergiant and main- 
sequence winds. The waist of this peanut is the Napoleon's 
Hat, which was previously thought to be an indepen- 
dent structure. The innermost circumstellar material 
is in the form of a cylindrical hourglass, 1 lt-yr in radius 
and 4 lt-yr lon g, connected to the pe anut by a thick 
equatorial disk (Sugerman et al. 2005a). With this 3D- 
model of the scattering dust the progenitor's mass-loss 
history can be reconstructed JSugerman et al.|2005b| . 
They found from the echo fluxes that from the interior 
hourglass to the bipolar lobes, the gas density drops 
from 1-3 to 0.03 cm -3 , while the maximum dust-grain 
size increases from 0.2 to 2 /im, and the silicate-to- 
carbonaceous dust ratio decreases, resulting in a total 
mass of ~1.7 Mq. The studies of |Sugerman et al.| 
( 2005a|b l are the most detailed three-dimensional cir- 
cumstellar dust reconstruction of any stellar object to 
date. 



4.4.2 SN 1987A light echoes scattered off 
interstellar dust 

It is of special interest to understand the structure 
of the ISM near SN 1987A in the LMC. This is a 
highly active starburst region with massive stars emit- 
ting strong stellar winds. The interaction of these 
high-energy outflows with the filamentary and spheri- 
cal structures of the ISM produce bubbles and super- 
bubbles in t he region which trigg er additional star for- 
mation (e.g., |Walborn et al.|1999t - LEs from SN 1987A 
allow the structure of these ISM supper-bubbles to be 
traced out in three dimensions. 

There are two nearly complete LE rings observed 
around SN 1987A, caused by two dust sheets roughly 
perpendicular to t he line-of-sight ~100 and ^400 pc 
in front of the SN (|Crotts|[T988l |Suntzeff et al.||1988b 



Gouiffes et al.||1988| |Couch et al.||1990[ ). |Xu eTa! 
( 1994 1 found two additional LEs at a much larger dis- 
tance of ~1000 pc in front of SN 198 7A. All of the 
above LEs were used by |Xu et al. ( 1995 1 to construct a 
detailed 3D-map of the ISM scattering dust. Compar- 
ing with Ha gas maps, Xu et al. ( 1995[ ) associate the 
ring at ~400 pc with reflection from the boundary of 
the super-bubble N157C. They found the curvature of 



this dust complex to coincide with LH 90, the source 
of the super-bubble. Additionally, the southeastern 
arc at a distance of ~1000 pc was found to align with 
a mas sive Ha filament discovered by |Meaburn et al.| 
(1984 1. They speculate that these two structures are 



the near and far sides of a giant super-bubble with a 
diameter of ~600 pc, which itself may have resulted 
from the merger of several smaller super-bubbles. 

4.4.3 Asymmetry of SN 1987A 

The two main SN 1987A LE ri ngs scattered off in- 
terstellar dust (see Section 4.4.21 provide an excellent 
opportunity to search for spectroscopic asymmetries 
in SN 1987A using the 3D spectroscopy technique de- 
tailed in Section [3.7| These LEs are very bright, allow- 
ing the observation of high signal-to-noise ratio spec- 
tra. A wealth of previous imaging of the LEs also ex- 
ists, making the observed apparent motion and hence 
dust inclination well known. Most importantly, how- 
ever, the original spectral and photometric history of 
SN 1987A are both well known, resulting in noambi- 
guity when modeling the LE spectra (Section 3.61. 

Sinnott et al. (in prep.) have obtained spectra of 
14 LEs of the inner echo ring of SN 1987A. In order 
to reduce the effects of nebular and stellar contamina- 
tion, they have also obtained follow-up sky-only obser- 
vations at the same position after the LE moved away 
using the same instrument configuration for the pur- 
pose of difference spectroscopy. In addition to provid- 
ing multiple viewing angles onto SN 1987A, the LEs 
which scatter off slightly different dust sheets while 
probing essentially the same viewing angle allows for 
a direct observational test of the LE spectra model- 
ing techniques described in Section |"3.6| After properly 
taking the dust filament extent and inclination into ac- 
count (see Figure [6]), they find evidence for asymmetry 
along a north-east to south-west axis, with little signs 
of asymmetry seen in the equatorial east-west direc- 
tions. The details of this asymmetry will appear in 
Sinnott et al. (in prep.). 

4.5 The 19th-century Great Erup- 
tion of rj Carinae 

r\ Carinae (n Car) is the most massive and most studied 
star in our Galaxy. It became the second-brightest star 
in the sky during its mid-19th century "Great Erup- 
tion," i n which it lost mo re than 10 Mq (Smith et al. 
2003b[ ). |Rest et al.| ( |2012[ ) discovered LEs of the Great 
Eruption, and subsequent spectroscopic follow-up re- 
vealed that its spectral type is most similar to those 
of G-type supergiants, rather than to F-type or earlier 
like typical luminous blue variable (LBV) outbursts. 
This raises doubts that traditional models involving 
opaque winds explaining LBV outbursts can fully ex- 
plain the Great Eruption. The absorption lines of the 
LE spectra indicate the emitting photosphere has an 
ejection velocity of ~ —200 km s _1 from a viewing an- 
gle perpend icular to the p rincipal axis of the Homuncu- 
lus Nebula ( Rest et al. 2012[). This is in agreement with 
velocities predicted by |Smith| ( |2006[ ). The n Car LE 
spectrum also shows a strong asymmetry in the Ca II 



www.pu blish . csiro. au/jo urnals /pasa 



13 



IR triplet, exten ding to a velocity of —850 km s 1 
( Rest et al.|2012 1. 

In recent years, some extragalactic non-SN tran- 
sients, (dubbed "SN impostors" ) , have been interpreted 



are much more likely caused by the intense and short 
burst of EUV-UV radiation associated with the shock 
breakout of the C as A SN itself ( |Dwek fe Arendt|2008| |. 
Kim et al. ( 2008 1 used the IR echoes to reconstruct the 



as analogues of the Great Eruption of r\ Car (e.g., Vink 3D-structure of the absorbing dust filaments 



20091, even though the Great Eruption is an extreme 



case in terms of energy, mass-loss, and duration. It is 
nevertheless surprising how different the LE spectra of 
the Great Eruption is compared to the SN impostors. 
Its spectral type is G2-G6, significantly later than all 
other SN impostors at peak, and it does not show any 
significant Hydrogen or Ca II IR triplet emission lines 
( |Rest et al.|2012[ ). 

4.6 Light Echoes of Ancient SNe in 
the LMC 

The first LEs of ancient or historic SNe were discovered 
in the LMC by |Rest et al.| ( |2005b| ). They found 3 
LE complexes associated with LMC SNRs 0509-675, 
0519-69.0, and N103B. Using measurements of the LE 
apparent motions, they determined the ages of these 
SNRs to be between 400 - 800 years. 

A spectrum of one echo associated with SNR 0509- 
675 reveals that the echo light is from the class of over- 
luminous Type la SNe (Rest et al. 2008a), the first 
time that an ancient SNe got typed based on its LE 
spectrum. This result is in excellent agreement with 
its classifi cation based on X-ray spectra of the SNR 
050 9-675 (|Hughes et aL|T995l |Badenes et al.|2008[ ) . 

|Hughes et al.| ( 1995 1 analyzed X-ray and radio data 
of the six SNRs in the LMC which have diameters 
less than 10 pc and are thus presumably the youngest 
in the LMC. They found three SNRs which have ele- 
ment abundances consistent with nucleosynthesis from 
Type II SNe, with the remaining three consistent with 
Type la nucleosynthesis. It should be noted that the 
three SNRs classified as Type la are also the SNRs 
with known LEs. This classification was confirmed for 
SNR 0509-675 ( |Rest et al.|2008a| ) and SNRs 0509-675 
and 0519-69.0 (Rest et al., in prep.). Only for the by 
far youngest core-collapse SN, SN 1987A, were LEs de- 
tected. Even though these are low-number statistics, 
the bias toward detecting LEs from Type la SNe can 
be explained by the fact that core-collapse SNe are 
typically much fainter than SN la. 

4.7 Cas A SN 

Cas A is the youngest (~330 years old) core-collapse 
SNR in our Galaxy (Stephenson & Green 2002 1. Krause 
et al. ( |2005 \ identified features in Spitzer Space Tele- 
scope IR images that changed with time. They iden- 
tified them as IR echoes, which are the result of dust 
absorbing the SN light, warming and re-radiating it 
at longer wavelengths. However, their main scientific 
conclusion that the cause of most if not all of these 
IR echoes was a series of recent X-ray outbursts from 
the compact object in the Cas A SNR turned out to 
be incorrect. They neglected to take into account that 
the apparent motion strongly depends on the inclina- 



The first scattered LEs of the Cas A SN were dis- 
covered by|Rest et al.| ( |2007||2008b| . At the same time, 
Krause et al. ( 2008a) spectroscopically followed up one 



of the IR echoes identified by the Spitzer Space Tele- 
scope, and found that the Cas A SN is most similar to 
the Type lib SN 1993J. This implies that the progen- 
itor of the Cas A SN was a red super giant that had 
lost most of its hyd rogen envelope before exploding 
( Krause et al.|2008a[ ) . 

Rest et al. ( 2011a[ ) obtained spectra of Cas A LEs 



from three different LEs spatially separated by degrees 
This means that each of the three spectra vie w th e 
Cas A SN from different viewpoints (see Section 3.71. 
They accounted for the effects of finite dust filament 
extent and inclination, and found that the He I A5876 
and Ha features of one LE are blue-shifted by an ad- 
ditional ~4000 kms -1 relative to the other two LE 
spectra and to the spectra of SN 1993J, indicating that 
Cas A was an intrinsically asymmetric SN. Data of the 
Cas A SNR in the X-ray and the optical also show that 
there is a dominant Fe-rich outflow in the same direc- 
tion (|Burrows et al.|20"05l|Wheeler et al.|2008[|DeLaney| 



et al. 12010 1, in excellent agreement with the LE data. 



This allows for the first time structure in the SNR to 
be directly associated with asymmetry observed in the 
explosion itself. 



4.8 Tycho's SN 



Tycho's SN of 1572, one of the last two naked eye 
SNe in the Galaxy, was classified as Type la based on 
its observed historical lightc ur ve and color ev o lution 

~ p006 l 



(Ruiz-Lapuente et al. 2004 1. Badenes et al 



finds that modeling the X-ray spectra of the Tycho 
SNR is most consistent with Tycho's SN being a nor- 
mal Type la, the first time that the subtype of Tycho 
was determined, even though only indirectly. 

Several LE groups of Tycho's SN were discovered 
by |Rest et al.] ( |2007| |2008b[ ), which opened t he door 
to spectroscopically classify this historical SN. |Krause| 
et al. | ( |2008b[ ) followed up one of these LEs, and classi- 
fied it as a normal SN la by comparing it qualitatively 
to a small number of SN LE templates constructed 
from spectropho tometry of nearby SNe , confirming the 
classification by Badenes et al. p006). It is notewor- 
thy that |Krause et al.| ( |2008b[ ) did not quantitatively 
compare the LE spectrum to a large set of comparison 
spectra templates, in particular they only compare the 
LE spectrum to a single underluminous SN la. |Rest| 
et al. ( 2008a I has shown that even comparing the high- 



tion o f the scattering d ust filament ( |Dwek fc Arendt 
2008] ) ( see also Section |3.3| . Instead, the IR echoes 



S/N LE templates to themselves only give conclusive 
classifications if a large number of templates are avail- 
able. 
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5 Conclusions 



Andrews, J. E., et al. 2011, ApJ, 731, 47 



In the last decade, the number of astrophysical ob- 
jects with observed LEs has dramatically increased. 
Newly-recognized LEs have been found around extra- 
galactic SNe as well as LEs from sources in the the 
Milky Way and Magellanic Clouds. In particular, LEs 
from nearby sources have proven to be fertile grounds 
for furthering our understanding of transients and vari- 
ables. The field has advanced from impressive geo- 
metric reconstructions of circumstellar and interstellar 
dust distributions, such as those produced from the 
study on SN1987A, to the realization of some of the 
early promise of improved geometric distances. The 
highly-unusual outbursting star V838 Mon had its dis- 
tance determined to 10% based on polarization mea- 
surements of the LEs by determining its geometric dis- 
tance with an accuracy of better than 10% l ] Sparks | 
et al. |2008" \ without any appeal to estimated luminosi- 
ties of possible counterparts. 

An exciting and promising new phase of SN sci- 
ence was ushered in with the discovery of whole LE 
systems around ancient /historic transients in the LMC 
and our Galaxy, including Cas A, Tycho, and r\ Car, 
where follow-up spectroscopy has allowed the direct 
comparison of the transient with its remnant. LEs 
have also allowed now been used, for the first time, to 
examine transients from the perspectives of their re- 
flection nebulae (scattering dust) and have provided 
observational spectroscopic constraints on the degree 
of asymmetry of outbursting objects. One as-yet un- 
realized opportunity provided by LEs is the ability to 
obtain spectroscopic time series of events and allow the 
reconstruction of the originally-unobserved lightcurve 
shapes of SNe. 

There is a great deal of discovery space left in LE 
research since only a small fraction of the sky has been 
searched for time- variable, low-surface brightness, non- 
stellar features. As currently envisaged, the next gen- 
eration of wide-field, time-d omain surveys like Pan - 
STARRS {Kaiser et afl2010"ll, PTF ( |Rau et al.|2009) , 
Skymapper ([Keller et al.||2007[ ), and ultimately LSST 
( Ivezic et al.|2008 l are being designed to study the ce- 
lestial sphere based dominantly on the variability of 
point-sources. A mode in which whole images of the 
sky are preserved with a cadence of weeks or months 
is likely to reveal an abundance of additional LE fea- 
tures and the history and perspectives they carry in 
their scattered light. 
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